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The  world’s  population  needs  new  sources  of  energy,  especially  those  that  are  clean  and  renewable. 
This  paper  provides  a  brief  introduction  to  wind  energy  and  the  types  of  existing  turbines,  which  are 
classified  using  the  orientation  of  the  rotation  axis.  Subsequently,  an  integral  analysis  is  performed  for 
vertical  axis  turbines.  The  known  variables  are  the  wind  speed,  the  type  of  blade,  the  radius  of  the  rotor 
and  the  angular  velocity.  The  fluid  velocity  and  the  angle  of  attack  on  the  blade  are  subsequently 
determined.  From  these  two  results,  the  lift  and  drag  forces  acting  on  the  blades  for  each  position  of 
the  rotor  are  calculated.  The  resultant  torque  and  power  generated  are  also  calculated  to  evaluate  the 
turbine  power  coefficient.  Due  to  the  rotation  and  the  robustness  of  this  type  of  turbine,  a  distortion 
in  the  flow  direction  occurs  in  its  vicinity.  The  flow  is  modeled  on  a  control  volume,  which  is  defined 
based  on  the  variation  in  the  wind  direction. 
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1.  Introduction 

Humanity  faces  a  major  challenge  in  meeting  energy  demands 
while  avoiding  damage  to  the  environment.  Energy  transported 
by  the  wind  is  part  of  the  solution  to  this  problem.  Wind  energy  is 
renewable  and  does  not  cause  any  harm  to  the  environment, 
thereby  enabling  the  generation  of  clean  electricity.  In  addition  to 
this  advantage,  the  cost  of  wind  energy  compared  to  the  cost  of 
conventional  systems  is  already  notably  competitive  [1],  Wind 
energy  produces  no  pollutants  in  contrast  to  a  thermal  power 
plant  that  produces  approximately  1  kg  of  carbon  dioxide  for  each 
kWh  produced. 

The  amount  of  energy  available  in  wind  varies  with  the 
seasons  and  the  time  of  day,  due  to  variations  in  wind  speed. 
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The  topography  and  roughness  of  the  soil  also  have  considerable 
influence  on  the  distribution  of  wind  speeds  in  one  given  place. 
Moreover,  the  amount  of  energy  extracted  from  the  wind  depends 
on  the  performance  characteristics,  time  of  operation  and  relative 
positions  of  the  wind  turbines. 

Wind  turbines  are  machines  used  to  extract  renewable  energy 
by  absorbing  part  of  the  wind’s  kinetic  energy  through  an 
aerodynamic  rotor.  The  wind  causes  the  turbine  blade  to  rotate 
around  its  axis,  converting  the  wind’s  kinetic  energy  into  mechan¬ 
ical  power,  which  is  subsequently  converted  into  electrical  power 
through  an  electrical  generator. 

The  amount  of  electricity  that  can  be  generated  using  the  wind 
depends  on  four  main  factors:  the  amount  of  air  passing  through  the 
turbine,  the  rotor  diameter,  the  type  of  blade  and  the  performance  of 
the  system. 

The  use  of  large  turbines  emerged  with  the  technology  of  the 
1970s  after  the  oil  crisis,  when  tax  incentives,  mostly  in  California 
(USA),  allowed  investments  in  these  systems.  Large  concentrations 
of  turbines,  called  Wind  Farms,  emerged.  After  the  removal  of  such 
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Fig.  1.  (a)  Horizontal  axis  wind  turbines  and  (b)  vertical  axis  wind  turbines. 

incentives,  the  pace  of  market  expansion  declined  but  resurged  at  the 
end  of  the  1990s.  In  this  new  millennium,  there  is  a  great  concern 
for  sufficient  energy  to  meet  the  increasing  consumption  of  both 
industrial  and  residential  users. 

A  wind  turbine  is  generally  classified  according  to  the  orienta¬ 
tion  of  its  rotor  axis  [2].  The  Horizontal  Axis  Wind  Turbine 
(HAWT),  shown  in  Fig.  la,  has  blades  that  rotate  in  a  plane 
perpendicular  to  the  main  direction  of  the  wind.  The  Vertical  Axis 
Wind  Turbine  (VAWT),  shown  in  Fig.  lb,  has  blades  that  rotate  in 
a  plane  parallel  to  the  direction  of  the  wind  [3], 

The  HAWT  is  the  most  commonly  used  wind  turbine  in  the 
market  to  produce  energy  on  a  large  scale,  but  it  requires  a  rotor 
axis  positioned  in  the  wind  direction  in  order  to  achieve  a  higher 
performance  efficiency,  especially  where  there  are  frequent  shifts 
in  the  wind  direction.  This  type  of  turbine  uses  a  tower  to 
maintain  the  turbine  components  at  an  optimum  height  where 
wind  speeds  are  higher.  However,  components  may  be  situated 
up  to  80  m  high,  creating  some  difficulty  when  performing 
inspections  and  maintenance. 

The  VAWT  has  the  ability  to  capture  wind  from  all  directions 
without  requiring  steering  systems  to  align  the  blades  with  the 
main  wind  flow.  Instead  of  a  tower,  it  is  also  possible  to  use 
support  cables  to  provide  structural  stability,  meaning  that  the 
elevation  of  the  rotor  can  be  lower  and  the  turbine  is  smaller  and 
less  expensive  to  construct.  With  the  blades  at  a  lower  elevation, 
the  generator  can  be  placed  at  the  ground  level,  thereby  facilitat¬ 
ing  inspection  and  maintenance  [4].  Conversely,  a  larger  base  area 
for  the  turbine  installation  is  required.  This  requirement  is  a 
major  disadvantage  in  agricultural  areas. 

The  first  VAWT  was  the  Panemonas  and  worked  with  the 
principle  of  drag  [2].  This  wind  turbine  probably  was  originated  in 
ancient  Persia,  China  or  Egypt.  However,  the  modern  VAWT, 
designed  using  modern  aerodynamic  concepts,  was  conceived 
by  the  Frenchman  D.G.M.  Darrieus  in  1920  [2,5].  There  are  several 
geometric  shapes  of  vertical  axis  wind  turbines,  which  can  be 
divided  into  three  basic  types,  called  Savonius  type,  Giromill  type 
and  Darrieus  type. 

The  operation  of  the  turbines  of  the  type  Savonius  occurs 
through  the  drag  force  on  its  blades,  which  are  two  cups  or  half 
drums  fixed  to  a  central  shaft  in  opposing  directions.  The  wind 
reaches  one  of  the  drums  and  the  drag  generated  on  this  causes 
the  axis  rotation.  Then  the  next  drum  reaches  the  position 
occupied  for  the  first,  and  the  drag  force  causes  the  axis  rotation 


a 


Fig.  2.  (a)  Savonius  type  VAWT,  (b)  Darrieus  type  VAWT,  and  (c)  Giromill 
type  VAWT. 

too.  This  process  continues  all  the  time  when  the  wind  blows. 
Fig.  2a  shows  a  Savonius  turbine. 

The  Darrieus  turbine  type  is  basically  lift  force  driven  wind 
turbine,  whose  long  and  flexible  blades  are  fixed  at  their  extremes 
and  deformed,  they  take  on  a  specific  configuration,  called  a 
troposkien,  which  minimizes  the  bending  stress  in  the  blades.  The 
Giromill  type  works  like  a  Darrieus  turbine,  but  uses  straight 
edges  and  constant,  became  the  case  to  be  produced  and  trans¬ 
ported,  reducing  costs.  Fig.  2b  and  c  shows  a  Darrieus  and  a 
Giromill  turbine  types,  respectively. 

A  disadvantage  of  VAWTs  is  the  fact  that  their  blades,  due  to  their 
rotation,  constantly  change  their  angles  of  attack  relative  to  the  wind 
direction,  resulting  in  alternating  and  time-dependent  forces.  This 
property  not  only  limits  their  efficiency  but  causes  structural 
vibrations.  Therefore,  the  choice  of  the  blade  geometry  is  of  para¬ 
mount  importance  when  improving  the  performance  of  the  turbines. 

According  to  the  literature  there  are  three  different  methodol¬ 
ogies  to  calculate  the  power  generated  by  a  vertical  axis  wind 
turbine  using  the  stream-tube  model:  the  single  stream-tube 
model,  the  multiple  stream-tube  model,  and  the  double-multiple 
stream-tube  model.  These  models  are  based  on  equating  the 
forces  on  the  rotor  blades  to  the  change  in  fluid  momentum  in 
the  flow  direction  through  the  rotor,  and  the  difference  between 
them  is  presented  as  follows. 

The  first  model  of  the  stream-tube  methodology  is  the  single 
stream-tube  model,  proposed  by  Templin  [6],  This  model  is  the 
simplest  model  because  this  assumes  a  constant  disk  velocity 
throughout  the  turbine.  In  other  words  the  speed  that  reaches  the 
blade  is  constant,  independent  of  the  angular  position  of  this 
blade,  because  this  model  predicts  a  uniform  flow  for  the  entire 
cross-section  (Fig.  3a).  The  forces  on  the  airfoil  blades  are  then 
computed,  using  this  uniform  velocity,  then  a  certain  error  is 
generated,  due  to  the  fact  that  the  speed  is  not  constant  in  each 
angular  position. 

A  new  model  was  created  in  order  to  consider  the  variation  in 
flow  cross-section  of  the  rotor,  multiple  stream-tube  model, 
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developed  by  Strickland  [7].  In  this  model  a  series  of  streamtubes 
are  assumed  to  pass  through  the  rotor  (Fig.  3b).  The  same  basic 
principles  which  were  applied  to  the  single  streamtube  are  applied 
to  each  of  the  multiple  streamtubes.  By  applying  the  momentum 
equation  is  each  streamtube  obtained  the  change  in  velocity  that 
reaches  the  blade  in  each  cross-section,  yields  a  more  realistic 
distribution  of  blade  forces,  so  this  model  provides  a  better 
representation  of  what  occurs  in  a  vertical  axis  wind  turbine. 

Paraschivoiu  [8]  proposed  a  more  sophisticated  model  called 
Double-Multiple  Stream-Tube  (Fig.  3c).  As  it  uses  two  actuator 
disks  placed  in  tandem  into  each  tube  of  the  multiple  array,  it  can 
also  predict  differences  between  the  upwind  and  downwind 
halves,  so  the  results  are  more  accurate. 

The  advancement  of  VAWTs  depends  on  the  success  of  current 
research  projects  and  the  performance  of  prototype  systems  under 
test.  Vertical  axis  turbines  are  the  focus  of  the  present  work. 

2.  Methodology 

A  physical  model  of  the  problem  is  first  presented  to  deter¬ 
mine  and  understand  the  variables  that  characterize  the  problem. 
Next,  a  mathematical  analysis  is  presented.  Finally,  simulations 
are  performed  for  various  values  of  the  parameters  in  order  to 
find  the  best  solution  to  this  problem,  i.e.,  maximum  efficiency. 

2.1.  Physical  model 

Fig.  4  presents  a  simplified  model  of  a  vertical  axis  turbine 
projected  on  a  horizontal  plane,  where  a  two-dimensional  analy¬ 
sis  is  performed. 

As  air  approaches  a  turbine,  it  is  disturbed  by  the  rotor, 
reducing  its  speed  from  the  free  stream  wind  speed  (Ua Q)  to  a 
lower  speed  U'  [9,10].  Moreover,  a  distortion  in  the  direction  of 
the  flow  occurs  in  the  vicinity  of  the  turbine,  mainly  due  to  the 
rotation  and  robustness  of  this  type  of  turbine.  In  other  words,  a 
streamline  that  tends  to  flow  in  one  direction  when  approaching  a 
turbine  tends  to  turn  in  the  same  plane  following  the  turbine. 


After  passing  through  the  turbine,  this  streamline  resumes  its 
original  direction.  Fig.  5  shows  a  magnified  representation  of  what 
occurs  near  the  turbine.  This  figure  represents  a  control  volume 
defined  by  planes  1  and  2  and  by  two  chosen  streamlines.  The  rotor 
is  immersed  in  this  control  volume.  A  second  control  volume  is 
delimited  by  planes  3  and  4  and  by  the  same  two  streamlines.  At 
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Fig.  3.  Different  types  of  stream-tube  models. 
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Fig.  6.  Resultant  velocity  ( V)  and  angle  of  attack  (a)  to  a  generic  position  6. 

plane  1  the  wind  speed  is  U (the  free  stream  speed),  and  the 
pressure  is  the  local  atmospheric  pressure.  Plane  2  is  considered  to 
be  sufficiently  distant  to  have  a  uniform  velocity  (U2)  and  a 
developed  flow  downwind  of  the  rotor,  where  the  pressure  is  also 
the  atmospheric  pressure.  Planes  3  and  4  define  the  inlet  and 
outlet  flows  over  the  turbine.  Due  to  the  speed  reduction,  it  is 
expected  that  the  pressure  in  plane  3  is  greater  than  atmospheric 
pressure,  and  the  pressure  in  plane  4  is  lower  than  atmospheric 
pressure,  due  to  the  viscous  effects  inside  the  rotor. 

A  generic  point  Q_  was  chosen  to  define  the  control  volumes 
used  in  the  analysis.  From  point  Q,  two  lines  were  drawn  tangent 
to  the  circle  defined  by  the  rotor,  thereby  defining  planes  3  and  4. 
Close  to  this  smaller  control  volume,  we  assume  that  the 
streamlines  take  the  shape  of  an  arc  over  the  rotor  with  their 
centers  at  the  point  Q,  For  example,  the  smaller  control  volume 
around  the  rotor  is  bounded  by  two  concentric  arcs  tangent  to  the 
rotor.  The  angle  y  is  formed  by  planes  1  and  3,  and  the  angle  f)  is 
formed  by  planes  3  and  4. 

Fig.  6  is  a  diagram  showing  the  wind  velocity  near  the  rotor  (LT), 
the  blade  velocity  ( coR )  and  the  resultant  velocity  (V).  The  figure 
shows  that  the  resultant  wind  velocity  for  any  angle  9  of  the  rotor  is 
the  sum  of  the  wind  speed  (U')  that  strikes  a  blade  and  the  relative 
velocity  between  the  fluid  and  the  rotating  blade  {coR),  where  co  is 
the  angular  speed  of  the  turbine  and  R  is  the  radius  of  the  rotor. 


Fig.  7.  (a)  Resultant  velocity  on  blade,  (b)  forces  acting  on  the  blade,  and 
(c)  resultant  force  in  the  direction  of  rotation. 


Knowing  the  resultant  velocity  of  the  fluid  over  the  blade,  the 
next  step  is  to  find  the  angle  of  attack  (a),  which  is  defined  as  the 
angle  between  the  blade  velocity  (ojR)  and  the  resultant  velocity 
( V ).  The  blade  velocity  ( coR )  is  placed  over  the  central  line  of  the 
blade  (called  the  chord  of  the  blade).  Fig.  6  shows  the  resultant 
velocity  and  angle  of  attack  (a)  for  a  generic  angle  6. 

Using  the  resultant  velocity  (Fig.  7a)  and  the  angle  of  attack, 
the  forces  acting  on  the  blade  are  determined.  These  are  the  drag 
force  (Fd)  and  the  lift  force  (F()  (Fig.  7b).  By  projecting  these 
forces  in  the  direction  of  the  blade  chord  (Fig.  7c),  the  tangential 
force  is  determined.  By  multiplying  the  tangential  force  by  the 
radius  R  of  the  turbine,  the  generated  torque  is  determined.  By 
multiplying  this  torque  by  the  angular  velocity,  the  generated 
power  is  obtained. 

To  finish  the  analysis,  the  remaining  variables  of  the  control 
volume  must  be  determined,  i.e„  the  U2  velocity,  pressures  in 
positions  3  and  4  (P3  and  P4,  respectively)  and  the  angles  ji  and  y. 

2.2.  Mathematical  model 

The  mathematical  model  for  the  described  problem  is  given  in 
the  following  steps,  which  are  defined  for  a  generic  position  9 
ranging  from  0  to  360 °,  thereby  giving  results  for  a  full  rotation  of 
the  rotor: 

1:  Estimate  the  wind  velocity  just  before  reaching  the  turbine 
(t/)  using  the  actual  wind  velocity  (Uoo) 

1/  =  ^  (1) 

2:  Calculate  the  resulting  fluid  velocity  on  the  blade  using 
cylindrical  coordinates.  First,  the  velocity  in  the  tangential 
direction  (l/9)  and  the  velocity  in  the  radial  direction  (Vr)  are 
calculated 

VB  =  coR+W  sin(0)  (2) 

Vr  =  U'  cos(9)  (3) 
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Finally,  the  resulting  velocity  is  calculated 

V=s/(Ve)2+(Vr)2  (4) 


3:  The  angle  of  attack  a  is  calculated  using  the  velocities  in  the 
radial  and  tangential  directions 


Vr 

a  =  arctan  — 
Vg 


(5) 


4:  The  Reynolds  number  Re  is  calculated  using  Eq.  (6)  [11] 


Re=^ 

P 


(6) 


where  p  is  the  fluid  density,  p  is  the  fluid  dynamic  viscosity  and 
c  is  the  blade  chord.  Because  the  Reynolds  number  and  the  angle 
of  attack  are  known  from  experimental  data  [12],  it  is  possible  to 
determine  the  drag  (Cd)  and  lift  (Q)  coefficients  for  the  blade  in 
any  position.  These  data  are  generated  as  follows:  for  different 
values  of  Reynolds  number,  Q  and  Cd  are  calculated  (from 
experimental  data)  for  angles  of  attack  ranging  from  0  to  180°. 
Next,  these  data  can  be  interpolated  in  order  to  determine  Cd  and 
Q  for  any  angle  of  attack  and  for  any  velocity. 

5:  Using  the  values  of  Q  and  Cd  for  the  resulting  velocity  and 
angle  of  attack,  the  coefficient  in  the  direction  of  the  axis  x' 
(C*')  is  calculated,  where  x'  is  the  direction  of  the  air  move¬ 
ment  (Fig.  4).  This  value  is  obtained  from  the  projection  of  Q 
and  Cd  on  x'  [13] 

Cxi  =(Cd  cos  a -Ci  sin  a)sin  6+(Cd  sin  a+C(  cos  a)cos  9  (7) 


6:  Using  Cx,  calculated  in  the  last  step,  the  velocity  of  the  wind 
that  reaches  the  turbine  is  estimated  and  compared  with  the 
velocity  estimated  in  step  1  to  verify  that  these  two  velocities 
have  the  same  value. 

Analyzing  the  control  volume  between  planes  1  and  2 
(Fig.  4),  we  can  say  that  the  speed  of  the  wind  at  the  turbine 
(this  step  is  called  U",  which  is  the  estimate  of  U'  in  each 
step)  is  the  average  between  the  input  velocity  (U^)  and  the 
output  velocity  ( U2 )  over  the  control  volume  [10] 

U"  =  U°°  +  U2  (8) 


The  force  in  the  x'  direction  is  given  by 

,  Cx,pU’2A 
F*  =  2 


(9) 


From  a  momentum  balance  [11],  the  force  in  the  x'  direction 
is  also  given  by 

Fx:  =  m(U00-U2)  (10) 


Isolating  U2  from  Eq.  (8)  and  using  Eqs.  (9)  and  (10),  we  have 


U"  = 


EL 


l+Cx-/4 


HD 


7:  The  error  between  U'  and  U"  is  given  by 
U'-U" 


(12) 


If  the  error  is  larger  than  a  given  residue  e,  one  must  return  to 
step  2  and  iterate  assuming  W  =  U"  [7], 

If  the  error  becomes  smaller  than  a  given  residue  e,  the  next 
step  is  taken. 

8:  Using  the  coefficients  (Q  and  Cd )  for  the  correct  speed,  the  lift 
force  (F|)  and  drag  force  (Fd)  (Fig.  6b)  are  calculated  [14] 


F,  =  C,\pAV2 


(13) 


Fd  =  Cd\pAV 2 


(14) 


9:  With  the  correct  values  of  F(  and  Fd,  the  resulting  forces  in  the 
normal  and  tangential  directions  are  obtained 

Fn  =  Fd  sin  a+F/ cos  a  (15) 


Fg  =  Fd  cos  a+Fi  sin  a 


(16) 


10:  The  next  step  is  to  calculate  the  generated  torque 

T=±FAe 

A-:  2  n 


(17) 


where  A 6  =  2n/N,  and  N  is  the  number  of  discrete  positions 
of  the  blade. 

1 1 :  The  generated  power  is  given  by 


P  =  Tco 


(18) 


12:  Using  the  generated  power,  the  power  coefficient  (Cp)  is 
calculated  [3,14-16] 

Cp=tA  09) 

1  max 

where  Pmax  is  the  maximum  power  that  can  be  obtained  from 
the  wind  [9,13] 

Pmax  =  JthU^ 

where  rh  =  pU^A  [9[.  Therefore 

Pmax=2P^^oo  (20) 

In  the  next  step,  the  remaining  variables  of  the  control 
volume  are  calculated. 

13:  The  resultant  force  on  the  system  in  the  smaller  control 
volume  (Fig.  4)  is  calculated  by  averaging  the  forces  acting  on 
the  turbine  for  each  complete  rotation.  These  forces  are  the 
force  in  the  direction  of  the  wind  motion,  FRmov,  and  the  force 
perpendicular  to  the  wind  motion,  FRperp 


r  Rmov 


_  -jL  Fg.  sindj+Fn,  COS0|- 


i  =  1 
N 


CRperp  —  ^2 
i  —  1 


2n 


Fg.  cosOi+Fn,.  sind,- 


2n 


Fr, 

and 

Fr  =  \/(FRm0V}  +(F/{p(,rp) 


Ad 


(21) 


(22) 


(23) 


Projecting  these  forces  in  the  x  and  y  directions,  we  obtain 
Frx  =  FRmov  cos (y+  0  +FRperp  sin^y  +  0  (24) 


FRy  =  FRmov  sin(j+  0  +FRperp  cos(j+  0 


(25) 


14:  With  this  resulting  force  in  the  wind  direction,  the  drag 
coefficient  over  the  rotor  can  be  defined  as 


Corotor  — 


PULA/2 


The  rotor  Reynolds  number  is  defined  as  [11] 


Re  rotor  — 


pujo 

P 


where  D  is  the  diameter  of  the  turbine. 


(26) 


(27) 
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15:  As  the  speed  of  the  wind  that  reaches  the  turbine  (U')  is  the 
averaged  input  speed  ( U ^  and  the  averaged  output  speed 
( U2 )  on  the  global  control  volume  [10],  U2  is  given  by 

U2  =  2U'-U00  (28) 

16:  Applying  the  Bernoulli  equation  [11]  between  positions 
1  and  3  (Fig.  4),  the  pressure  P3  is  obtained 

P3  =  Pa,m+^(Ul-U'2)  (29) 

where  PaCm  is  the  atmospheric  pressure. 

Likewise,  using  the  Bernoulli  equation  between  positions  4 
and  2  (Fig.  4),  the  pressure  P4  is  also  obtained 

P4  =  Patm+^(U22-U'2)  (30) 

17:  From  a  linear  momentum  balance  [11  j  applied  on  the  smaller 
control  volume  (between  planes  3  and  4)  and  projecting  the 
forces  in  the  x  and  y  directions,  we  obtain  two  equations 
involving  two  unknown  variables,  p  and  y 

P  r +P3  A  3  +P4  A  4  =  0  (31) 

Projecting  the  forces  in  the  x  and  y  directions,  we  obtain 
Frx  =P4  sin(90D-y-/i)A-P3  cos(y)A  (32) 


FRy  =  P4  cos(90 °-y-fi)A-P3  sin(y)A  (33) 

The  variables  p  and  y  can  be  obtained  from  Eqs.  (32)  and  (33), 
as  all  other  variables  have  already  been  determined. 
When  we  equate  Eq.  (24)  to  Eq.  (32)  and  Eq.  (25)  to 
Eq.  (33),  we  obtain  a  nonlinear  system  of  equations.  To  solve 
this  system,  we  can  use  an  iterative  method.  First,  we 
estimate  an  initial  value  /?'  for  p.  Next,  we  calculate  p  as  a 
function  of  p' 


P  =  2  arccos 


F2r+(P4A)2+(P3A)2  \ 

2piA{F«-Pi^)+2p4AF«J 


(34) 


Eq.  (34)  is  solved  until  the  angle  p  converges  assuming  /('  =  p. 
After  calculating  /?,  we  calculate  y  directly 


y  =  arctan 


f  Fr  cos(P/2)+P4A  cos(/?/2)— P3A‘ 

V  /V^irt^FP^WZ)  . 


(35) 


Certain  variables,  which  are  the  resultant  torque  (T),  the 
power  coefficient  (Cp),  the  average  velocity  at  the  turbine 
((/)  and  the  velocity  at  position  2  (l/2),  the  pressure  acting 
upward  (P3)  and  downward  (P4)  on  the  turbine  and  the 
angles  y  and  p  due  to  the  distortion  of  the  streamlines 
around  the  turbine,  will  be  plotted  as  a  function  of  tsr 
(tip  speed  ratio),  a  dimensional  variable  which  is  the  ratio 
of  the  rotational  speed  of  the  blade  to  the  free  wind  speed 
[15,16] 


tsr  = 


coR 


(36) 


With  this  variable,  we  can  calculate  the  forces,  torque  and  power 
for  various  blade  rotational  speeds  and  wind  speeds  to  deter¬ 
mine  the  point  of  the  maximum  efficiency. 

It  is  important  to  stress  the  viscous  effects  are  taken  into  account 
in  the  integral  formulation  as  well  as  the  experimental  results 
are  used  together  the  present  theoretical  model. 


3.  Results  and  discussion 

The  following  data  were  assumed  in  order  to  perform  the 
simulations  and  to  obtain  the  results  presented  below.  These 


values  were  chosen  arbitrarily  to  illustrate  the  mathematical 
methodology  presented  (Table  1 ). 

Figs.  8  and  9  show  the  lift  coefficient  and  the  drag  coefficient 
as  a  function  of  the  angle  of  attack  of  the  chosen  profile,  NACA 
0012.  This  coefficients  was  obtained  experimentally  [12]. 

The  following  figures  present  the  solutions  for  the  variables  of 
the  problem  using  the  presented  mathematical  model.  Fig.  10 
shows  the  velocity  U'  at  which  the  wind  reaches  the  turbine.  This 
velocity  was  obtained  by  applying  the  iterative  method  (steps  2-7 
of  the  mathematical  model). 

Fig.  10  shows  that  the  velocity  U1  of  the  wind  that  reaches  the 
blade  follows  a  periodic  function.  When  the  chord  of  the  blade 
is  in  the  direction  of  the  air  flow  (90°  and  270°),  the  velocity 
approaches  the  free  wind  velocity  (10  m/s).  When  the  chord  of 
the  blade  is  perpendicular  to  the  blade  motion  (0°  and  180  ),  the 
velocity  of  the  air  at  the  turbine  is  minimal.  This  finding  indicates 
that  as  the  blade  chord  approaches  this  position,  the  air  velocity 
at  the  blade  decreases.  It  is  also  shown  that  the  maximum  value 
of  the  air  velocity  is  the  same  for  all  values  of  tsr,  and  as  the  tsr 
increases  (increasing  the  angular  velocity),  the  minimum  value 
for  U'  also  increases,  i.e.,  the  deceleration  of  the  wind  is  reduced. 

The  resulting  velocity  on  the  blade  and  angle  of  attack  are 
calculated  and  are  shown  in  Figs.  11  and  12.  The  velocity  on  the 
blade  reaches  its  maximum  value  at  9  =  90°  and  its  minimum 
value  at  9  =  270  ,  when  the  wind  and  the  blade  move  in  the  same 
and  opposite  directions,  respectively. 

Increasing  the  tsr  (thus  increasing  the  rotor  angular  velocity) 
increases  the  resulting  velocity  of  the  turbine  in  the  same  manner 
for  all  positions,  as  indicated  by  a  displacement  of  the  curve. 

Note  that  in  Fig.  12,  as  the  rotation  of  the  rotor  increases,  the 
angle  of  attack  decreases  [15].  This  occurs  because  the  fluid 
velocity  in  the  tangential  direction  also  increases.  Note  also  that 


Table  1 

Values  of  variables  used  for  simulation. 


Variable 

Value 

Free  wind  speed 

UOQ  =  10  m/s 

Turbine  radius 

R=0.64  m 

Chord  length  (Profile  NACA  0012) 

c=3.2  cm 

Tip  speed  ratio 

From  4  to  1 8 

Position  of  the  blade 

From  0  to  360° 

Density 

p  =  1.19  kg/m3 

Viscosity 

p  =  1.84F-5  kg/(m  s) 

Atmospheric  pressure 

Patm  =  1.01325E5  Pa 

Fig.  8.  Experimental  lift  coefficient  (C/)  as  a  function  of  the  angle  of  attack  (a)  of 
the  profile  NACA  0012  for  various  Reynolds  number  [12]. 


R.  Romao  da  Silva  Melo,  A  da  Silveira  Neto  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  4809-4817 


4815 


Fig.  9.  Experimental  drag  coefficient  (Cd)  as  a  function  of  the  angle  of  attack  (a)  of 
the  profile  NACA  0012  for  various  Reynolds  number  [12]. 


Fig.  10.  Wind  velocity  at  the  turbine  (17')  for  various  tsr  values  as  a  function  of  0. 


the  angle  of  attack  is  zero  when  the  blade  is  in  the  direction  of  the 
wind  velocity,  as  expected. 

Fig.  13  shows  the  resultant  force  in  the  tangential  direction. 
Note  that  the  tangential  force  is  at  a  maximum  when  the  chord 
blade  is  close  to  the  perpendicular  direction  of  the  wind  (near 
0  =  0  and  180°);  therefore,  in  these  positions,  the  torque  is  also  at 
a  maximum.  When  the  chord  blade  is  close  to  the  normal 
direction  of  the  wind  (0  =  90  and  270°),  the  tangential  force  is 
minimal  and  even  negative  at  some  positions.  However,  it  is 
observed  that  during  each  cycle,  the  tangential  force  is  generally 
positive,  and  the  total  torque  is  therefore  also  positive. 


Fig.  14  shows  the  resultant  normal  force.  This  figure  shows 
that  the  absolute  value  of  the  normal  force  is  larger  when  the 
chord  blade  is  close  to  the  perpendicular  direction  of  the  wind 
(near  9  =  0  and  180°),  and  smaller  when  the  chord  blade  is  close 
to  the  parallel  direction  of  the  wind  (0  =  90°  and  270°).  As 
the  angular  velocity  increases,  the  normal  force  also  increases. 
The  normal  force  does  not  influence  the  resultant  torque  but  is 
important  in  determining  the  mechanical  efforts  on  the  rotor. 

Fig.  15  shows  the  torque  (T)  generated  by  the  turbine,  which 
was  obtained  using  the  average  of  the  tangential  force  acting  on  it 
(Eq.  (17)).  The  torque  is  bigger  for  lower  values  of  tsr,  due  to  the 
fact  that  tangential  forces  also  are  bigger  for  this  tsr  values. 

We  can  also  obtain  the  power  (Eq.  (18))  and  the  power 
coefficient  Cp  (Eq.  (19)),  which  is  shown  in  Fig.  16. 
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Fig.  15.  Torque  on  the  turbine  as  a  function  of  tsr. 


tsr 


Fig.  16.  Power  coefficient  of  the  turbine  as  a  function  of  tsr. 

The  Cp  curve  obtained  obeys  the  characteristics  of  this  type  of 
turbine.  As  the  tsr  increases,  the  Cp  increases,  reaching  a  maximum 
before  decreasing.  In  this  simulation,  the  maximum  Cp  is  near  0.35 
for  a  tsr  near  11 .  It  can  also  be  observed  that  the  power  coefficient 
does  not  reach  the  Betz  limit,  which  is  0.59  [10],  as  expected. 

Fig.  17  shows  a  comparison  between  the  three  velocities:  the 
free  stream  wind  velocity  ( U «,),  the  average  velocity  at  the  turbine 
(Eq.  (11))  and  the  velocity  at  position  2  (Fig.  5)  of  the  larger  control 
volume  (U2,  Eq.  (28)). 

The  velocity  of  the  wind  reaching  the  turbine  increases  as 
the  turbine  rotation  increases,  and  when  the  angular  velocity 
approaches  infinity,  the  average  velocity  approaches  the  free 
stream  wind  velocity.  The  velocity  at  point  2  is  less  than  V  and 
is  directly  proportional  to  the  velocity  at  the  turbine.  The  curve 
shows  that  the  air  velocity  at  the  turbine  is  the  average  of  the  free 
stream  wind  velocity  and  the  velocity  at  position  2. 

Fig.  18  shows  the  averaged  forces  acting  upon  the  turbine 
during  one  complete  rotation.  After  averaging  the  forces  during 
one  complete  rotation,  we  see  that  the  force  in  the  direction 
perpendicular  to  the  motion  of  the  fluid  is  zero,  leaving  only  the 
force  in  the  direction  of  the  air  movement.  This  force  increases 
linearly  with  the  angular  velocity. 


10.5 


10 


8 


7.5  - ' - ' - ' - 

6  9  12  15  18 

tsr 

Fig.  17.  Resulting  velocities  as  a  function  of  various  tsr  values. 


Fig.  18.  Resulting  average  forces  for  various  tsr  values. 


Fig.  19.  Drag  coefficient  over  the  rotor  as  a  function  of  the  Reynolds  number. 

Fig.  19  shows  the  drag  coefficient  over  the  rotor  as  a  function 
of  the  Reynolds  number.  The  drag  coefficient  decays  nonlinearly, 
as  expected.  This  decay  occurs  because  the  drag  coefficient  on  the 
rotor  is  proportional  to  the  drag  force,  which  increases  linearly 
and  is  inversely  proportional  to  the  square  of  the  wind  speed. 

Fig.  20  shows  the  pressure  acting  upward  (P3)  and  downward 
(P4)  on  the  turbine  and  the  atmospheric  pressure.  It  is  observed 
that  the  pressure  P3  before  the  turbine  is  above  atmospheric 
pressure  and  the  pressure  P4  after  the  turbine  is  below  atmo¬ 
spheric  pressure.  This  is  due  to  the  drag  generated  by  the  turbine. 
Note  also  that  the  pressure  decays  before  the  turbine  and 
increases  after  the  turbine  as  the  tsr  increases.  When  the  speed 
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Fig.  20.  Pressures  distributions  for  various  tsr  values. 


Fig.  21.  Pressures  distributions  for  various  tsr  values. 


Fig.  22.  Pressures  distributions  for  various  tsr  values. 


ratio  approaches  infinity,  the  pressure  before  and  after  the  turbine 
approaches  the  atmospheric  pressure. 

Figs.  21  and  22  show  the  angles  (5  and  y  due  to  the  distortion  of 
the  streamlines  around  the  turbine.  It  is  expected  that  as  the 
angular  velocity  decreases,  the  streamlines  display  less  distortion, 
and  when  the  rotation  angle  is  zero,  both  angles  approach  zero. 


Otherwise,  the  higher  the  angular  velocity,  the  greater  the 
distortion  of  the  streamlines. 


4.  Conclusion 

Wind  energy  is  part  of  the  solution  to  the  world’s  energy 
demands.  Computational  analysis  for  wind  turbines  is  notably 
versatile,  as  it  avoids  the  costs  of  experiments  and  projects  in 
physical  facilities,  and  it  can  easily  simulate  and  handle  variable 
changes  to  seek  the  best  solutions. 

With  a  full  analysis  for  given  wind  conditions,  it  is  possible  to 
choose  a  good  aerodynamic  profile  for  a  turbine.  The  angular 
velocity  of  the  axis  can  be  chosen  such  that  the  turbine  operates 
near  its  maximum  efficiency. 

With  the  coefficient  of  total  drag  on  the  rotor  and  the  angles  of 
distortion,  it  is  possible  to  perform  a  global  analysis  of  the  flow  in 
a  wind  farm.  We  can  determine  the  minimum  distance  required 
between  two  turbines  to  guarantee  that  the  air  flow  passing 
through  the  first  turbine  approaches  the  velocity  and  pressure  of 
the  free  wind  before  passing  through  the  next  turbine. 

Nevertheless,  the  integral  analysis  requires  empirical  inputs 
for  each  type  of  rotor  blade.  In  order  to  avoid  this  requirement, 
computational  fluid  dynamics  methods  can  be  useful.  This  subject 
is  the  next  step  and  the  focus  of  the  continuation  of  the  present 
research. 
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